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Abstract 

The  settling  velocities  of  undisturbed  macroscopic  aggregates  known  as  marine  snow  were  mea¬ 
sured  with  SCUBA  in  surface  waters  off  southern  California  and  analyzed  as  a  function  of  aggregate 
size,  mass,  and  density*  The  mean  ^ttling  velocity  was  74±39md“'for  aggregates  ranging  from 
2.4  to  75  mm  in  maximum  lengtfirSinking  rates  in  the  field  varied  exponentially  with  aggregate 
size  and  dry  weight  and  were  consistently  up  to  four  times  slower  than  rates  measured  in  the 
IgboratoQj  ^ "  '  -'ur,  ' 

^^he  excess  densities  of  the  80  aggregates  examined  were  calculated  from  volume  and  dry  weight 
and  ranged  over  four  orders  of  magnitude  with  a  median  oi[X4  x  10"^  g  cm"’'.  Aggregates  of 
marine  snow  sank  more  slowly  than  predicted  for  either  solid  or  perous  spheres  of  equivalent 
volume  and  density,  although  their  velocities  were  within  the  range  expected  for  equivalent  sinking 
prolate  ellipsoids.  No  relationships  between  settling  velocity  and  either  excess  density  or  particle 
shape  were  found.  Drag  coefficients  of  marine  snow  were  also  higher  than  predicted  by  theory  for 
spheres  of  equivalent  volume  and  density.  These  deviations  from  theoretical  expectations  may  be 
partially  explained  by  errors  in  the  estimation  of  the  excess  densities  of  aggregates.  Variability  in 
the  densities  of  the  heterogeneous  primary  particles  comprising  marine  snow  (fecal  pellets,  clay- 
mineral  particles,  phytoplankton,  molts,  etc.)  and  the  potential  for  buoyancy  regulation  by  indi¬ 
vidual  phytoplankton  cells  inhabiting  aggregates  make  determination  of  excess  density  especially 
problematic,  .  'i  >  .  t 

Large,  rapidly  sinking  particles  are  far  less  (Silwr  and  Alldredge  1981;  Taguchi  1982; 

numerous  than  fine  suspended  particles  in  Gorsky  et  al.  1983)  and  the  field  (Shanks 

the  pelagic  zone  of  the  sea.  However,  recent  and  Trent  1 980)  and  estimated  from  in  situ 

studies  suggest  that  large  particles  are  pri-  photographs  (Billett  et  al.  1983;  Lampitt 

rnarily  responsible  for  the  vertical  transport  1985;  Asper  1987),  the  relationships  be- 

of  biogenic  materials  and  eletnents  through  tween  settling  velocity  and  aggregate  prop- 

the  water  column  (see  Fowler  and  Knauer  erties  such  as  size,  mass,  and  density  are  not 

1986).  Although  fecal  pellets  are  a  major  known.  Accurate  estimates  of  flux  based  on 

source  of  large,  rapidly  sinking  particles  in  particle  abundances  in  nature  require  in- 

Svxme  regions  of  the  world’s  oceans  (Angel  formation  on  size-specific  sinking  rates. 

1984),  in  most  other  areas  flocculent  amor-  Moreover,  data  on  the  settling  behavior  of 

phous  aggregates  >0.5  mm  in  diameter,  aggregates  are  necessary  in  order  to  under¬ 
known  as  marine  snow,  appear  to  be  the  stand  the  quantitative  impact  of  processes 

major  source  ofparticulate  flux  (Fowler  and  that  alter  aggregate  size  and  density  (con- 

Knauer  1986).  sumption  by  grazers,  dissolution,  disaggre- 

Sinking  aggregates  transfer  energy  and  gation,  etc.)  on  the  biogeochemical  pro¬ 
nutrients  from  the  surface  mixed  layer  into  cesses  associated  with  particles  in  the  water 

the  ocean  interior  and  to  the  seafloor.  The  column.  In  this  paper  we  present  the  first 

rate  at  which  this  transfer  occurs  depends  data  on  the  settling  behavior  of  undisturbed 

primarily  on  the  composition,  size,  abun-  aggregates  of  marine  snow  in  situ  as  a  func- 

dance,  and  sinking  rate  of  the  particles.  Al-  tion  of  aggregate  size,  mass,  and  density.  We 

though  settling  velocities  of  marine  snow,  discuss  variations  in  this  settling  behavior 

ranging  from  1  to  370  m  d""',  have  been  from  predicted  theory  and  provide  evidence 

measured  directly  in  both  the  laboratory  that  laboratory  studies  substantially  over¬ 

estimate  sinking  rates  of  marine  snow  in 
nature. 
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Methods 

Theoretical  considerations— a  par¬ 
ticle  settles  at  a  constant,  or  terminal,  ve¬ 
locity  through  a  fluid,  the  force  of  gravity 
pulling  the  particle  downward  is  exactly  bal¬ 
anced  by  the  retarding  drag  force  of  the  fluid 
flowing  around  the  particle.  The  resulting 
force  balance  equation  takes  the  form 

v{p,-  p;)g=yiCoApfiP  (1) 

where  V  is  the  volume  of  the  particle  in  cm^. 
Pa  the  density  of  the  particle  in  g  cm  “^,  pr 
the  density  of  the  fluid  in  g  cm“^,  g  the 
acceleration  due  to  gravity  or  980  cm  s' 

A  the  maximum  cross-sectional  area  of  the 
particle  perpendicular  to  the  direction  of 
sinking  in  cm^,  U  the  settling  velocity  of  the 
particle  in  cm  s”‘,  and  Cp  (the  drag  coeffi¬ 
cient)  a  dimensionless  number.  For  spher¬ 
ical  particles,  Co  is  a  function  of  Reynolds 
number  (Re  =  dCIv,  where  d  is  particle  di¬ 
ameter  in  cm  and  v  is  the  kinematic  vis¬ 
cosity  of  the  fluid).  At  Reynolds  numbers 
<0.5,  Stokes’  expression  for  a  settling  sphere 
is  valid  and  s  24/Re  (White  1974).  For 
nonspherical  particles,  such  as  marine  snow, 
however,  the  drag  coefficient  becomes  a 
complicated  function  of  both  shape  and 
Reynolds  number.  This  function  must  usu¬ 
ally  be  determined  empirically  (see  Tietjens 
1957;  Graf  1971;  Komar  and  Reimers 
1978).  Rearranging  Eq.  1  yields 

U  =  (2g^pV/p,CoAy’  (2) 

where  Ap  =  (pa  ~  p/).  From  Eq.  2  we  would 
expect  the  terminal  settling  velocity  of  an 
aggregate  of  marine  snow  to  be  a  function 
of  the  excess  density  of  the  particle,  the  ratio 
of  volume  and  projected  area  (in  the  direc¬ 
tion  of  settling),  and  the  drag  coefficient.  We 
measured  the  physical  properties  of  marine 
snow  appropriate  for  elucidating  these  re¬ 
lationships. 

Since  the  drag  coefficient  and,  thus,  set¬ 
tling  velocity  are  a  function  of  shape,  we 
also  required  a  coefficient  that  adequately 
describes  the  shape  of  an  aggregate  and  how 
it  departs  from  spherical.  Most  previous 


studies  of  the  effects  of  shape  on  the  settling 
of  natural  particles,  including  sediment 
grains  (Komar  and  Reimers  1978)  and  small 
animals  such  as  foraminifera  (Fok-Pun  and 
Komar  1983),  have  used  the  Corey  shape 
factor  (CSF)  (Albertson  1 953;  see  McNown 
and  Malaika  1950): 


CSF 


Ds 

(D.Dy^ 


(3) 


where  D^,  D,,  and  £>,  are  the  smallest,  largest, 
and  intermediate  axial  diameters  of  the  par¬ 
ticle.  For  a  spherical  shape,  CSF  =  1 .  The 
further  from  spherical  a  particle  becomes, 
the  closer  to  0  CSF  becomes.  McNown  and 
Malaika  (1950)  found  that  the  ratio  of  the 
principal-axis  lengths  was  the  best  predictor 
for  determining  the  effect  of  shape  on  the 
settling  velocities  of  small  ellipsoids,  cones, 
cylinders,  and  prisms  up  to  Re  of  ~  10. 

The  simplest  way  to  consider  the  effects 
of  shape  on  settling  velocity  is  to  consider 
a  sphere  deformed  to  form  a  new  shape  of 
equivalent  volume  and  density,  encounter¬ 
ing  a  different  resistance  from  the  fluid  flow 
past  it.  This  resistance  differs  from  that  of 
the  original  sphere  by  k,  the  coefficient  of 
form  resistance  (also  called  the  dynamic 
shape  factor)  where 


k  = 


U 


(4) 


with  (7j  the  sinking  velocity  of  a  sphere  of 
equivalent  volume  and  density  and  U  the 
settling  velocity  of  the  particle.  Theoretical 
derivation  of  k  and  further  discussion  of  the 
effects  of  particle  shape  on  it  can  be  found 
elsewhere  (McNown  and  Malaika  1950; 
Davis  1979;  Hutchinson  1967).  For  ellip¬ 
soids  of  regular  shape,  k  varies  as  a  direct 
function  of  CSF.  Most  prolate  (largest  axial 
diameter  parallel  to  ffow)  ellipsoids  and  all 
oblate  (largest  axial  diameters  perpendicu¬ 
lar  to  flow)  ellipsoids  have  a  k  >  1  and  sink 
more  slowly  than  equivalent  spheres.  The 
further  from  spherical  the  particles  become 
the  larger  k  becomes  (McNown  and  Malaika 
1950). 

Empirical  measwrcmcnts— Settling  ve¬ 
locities  of  individual  aggregates  of  marine 
snow  were  measured  both  in  situ  and  in  the 
laboratory.  Field  measurements  were  made 
at  depths  of  10-15  m  during  25  SCUBA 
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dives  in  the  San  Pedro  Basin,  California  (4 
km  due  east  of  Two  Harbors,  Catalina  Is¬ 
land),  in  March  and  June  1986  and  during 
two  dives  in  the  Santa  Barbara  Channel  (6 
km  due  south  of  the  University  of  Califor¬ 
nia,  Santa  Barbara)  during  July  1986.  All 
dives  were  in  late  morning.  We  measured 
the  settling  velocities  of  80  undisturbed  ag¬ 
gregates  of  marine  snow  relative  to  a  spot 
of  dilute,  neutrally  buoyant  fluorescein  dye 
placed  3  cm  below  each  particle.  The  pow¬ 
dered  dye  was  mixed  with  seawater  at  depth 
so  as  to  produce  a  solution  of  the  same  den¬ 
sity  as  the  seawater  surrounding  the  aggre¬ 
gates.  A  SCUBA  diver  delivered  the  spot  of 
dye  (a  few  millimeters  in  diameter)  from  a 
dye  ejector  consisting  of  a  hypodermic 
needle  mounted  on  a  1-m-long,  2.4-cm-di- 
ameter  plastic  pipe  (Fig.  1).  The  stock  of 
this  dye  ejector  rested  against  the  diver’s 
shoulder  to  stabilize  the  ejectf^r  and  reduce 
vibration  at  its  tip. 

The  3-cm  distance  was  determined  by  sit¬ 
ing  the  bottom  of  the  aggregate  at  the  tip  of 
a  thin  wire  attached  to  the  ejector  exactly  3 
cm  above  the  release  point  of  the  dye  (Fig. 
la).  The  dye  spot  congealed  and  stabilized 
in  shape  and  size  within  1-2  s  after  delivery. 
Several  minutes  were  required  for  it  to  dif¬ 
fuse.  With  care,  a  diver  could  deliver  a  spot 
of  dye,  slowly  pull  the  ejector  back,  and  hang 
suspended  2-3  m  from  the  aggregate  with¬ 
out  producing  turbulence  near  the  aggre¬ 
gate.  Turbulence  was  easily  detected  by 
observing  deformation  of  the  dye  spot.  A 
second  diver,  located  7-8  m  away,  mea¬ 
sured  with  a  stopwatch  the  time  required 
for  the  particle  to  sink  to  the  dye  spot,  re¬ 
lying  on  visual  cues  from  the  first  diver.  This 
second  diver  then  photographed  the  still  un¬ 
disturbed  aggregate  to  provide  a  record  of 
its  size  and  shape  using  a  Nikonos  camera 
with  a  1 :  1  or  3 :  1  close-up  extension  tube 
and  framer.  The  aggregate  was  then  col¬ 
lected  in  a  6-ml  plastic  cylinder  for  later 
determination  of  aggregate  dry  weight. 

A  single  photograph  was  taken  in  a  plane 
parallel  to  the  direction  of  sinking  (i.e.  a  side 
view)  by  orienting  the  particle  at  a  prede¬ 
termined  focal  distance  from  the  camera  lens 
demarcated  by  a  thin  wire  frame.  A  trans¬ 
parent  metric  ruler  photographed  in  the 
frame  yielded  scale.  Turbulence  produced 


by  moving  the  camera  after  the  first  pho¬ 
tograph  altered  the  aggregate’s  shape  and  in 
some  cases  threatened  to  break  the  particle 
apart.  These  factors  thwarted  our  attempts 
to  obtain  additional  photographs  in  other 
planes.  Direct  observation  of  many  aggre¬ 
gates  in  the  field  indicated,  however,  that 
most  sinking  aggregates  are  shaped  like 
comets  or  prolate  ellipsoids  and  are  nearly 
symmetrical  about  their  polar  axes  (parallel 
to  the  direction  of  sinking).  Thus,  a  pho¬ 
tograph  oriented  parallel  to  the  plane  of 
sinking  allowed  for  the  most  accurate  esti¬ 
mation  of  volume. 

The  sc'tling  velocities  of  37  different  ag¬ 
gregates  were  determined  in  the  laboratory 
for  comparison  with  field  measurements. 
Aggregates  of  marine  snow  collected  from 
lO-m  depths  in  the  Santa  Barbara  Channel 
were  allowed  to  sink  slowly  out  of  the  open 
bottom  of  their  collecting  cylinders  into  the 
top  of  an  8-cm-diameter,  2,U00-ml  gradu¬ 
ated  cylinder  containing  seawater  from  the 
depth  of  collection.  Both  the  samples  and 
the  graduated  cylinder  were  maintained  in 
an  environmental  chamber  at  ambient  sea¬ 
water  temperature  ( 1 6°C),  and  settling  rates 
were  determined  within  4  h  of  sample  col¬ 
lection.  Periodic  injection  of  fluorescein  dye 
into  the  settling  chamber  indicated  that  no 
convective  currents  were  present.  We  timed 
the  sinking  of  each  aggregate  past  external 
markings  on  the  graduated  cylinder  with  a 
stopwatch.  Each  aggregate  was  then  re¬ 
moved  from  the  chamber  with  a  pipette  and 
placed  on  a  filter  for  dry  weight  determi¬ 
nation.  Dry  weight  was  used  as  the  property 
of  marine  snow  potentially  least  affected  by 
collection  and  handling. 

Physical  characteristics  of  each  sinking 
aggregate  were  determined  as  follows— 

Dry  weight:  Each  individual  aggregate  was 
filtered  with  3  ml  of  surrounding  seawater 
onto  a  preweighed  Nuclepore  filter  (2.4-cm 
diameter,  0.40-/xm  pore  size),  rinsed  quickly 
with  distilled  water  and  dried  for  a  mini¬ 
mum  of  24  h  in  a  desiccator.  Filters  were 
reweighed  to  the  nearest  0. 1  ^g  on  a  Cahn 
electrobalance  and  the  dry  weights  were  cor¬ 
rected  with  3-ml  seawater  blanks.  Reweigh¬ 
ing  of  some  filters  10  d  after  placement  in 
the  desiccator  indicated  that  all  water  was 
lost  within  the  first  24  h. 

I  “  ‘  '  ■  -4 
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Diameter;  The  maximum  diameter  per¬ 
pendicular  to  the  direction  of  sinking  was 
measured  to  the  nearest  0. 1  mm  on  the  pho¬ 
tographs. 

Maximum  length:  The  maximum  length 
was  measured  to  the  nearest  0. 1  mm  on  the 
photographs. 

Projected  area:  Drag  on  a  sinking  particle 
may  be  correlated  with  the  maximum  cross- 
sectional  area  perpendicular  to  the  direction 
of  sinking.  We  calculated  this  projected  area 
with  the  measured  maximum  cross-sec¬ 
tioned  diameter  assuming  that  this  area  was 
a  circle.  Photographs  taken  perpendicular 
to  the  direction  of  sinking  would  have  yield¬ 
ed  the  projected  area.  Photographing  small, 
fragile  objects  from  directly  above  or  below, 
however,  is  nearly  impossible  for  a  diver 
floating  in  open  water.  Moreover,  such  pho¬ 
tographs  would  have  yielded  underesti¬ 
mates  of  volume  since  most  aggregates  were 
elongated  along  their  polar  axes. 

We  estimated  the  error  involved  in  cal¬ 
culating  the  projected  area  from  the  maxi¬ 
mum  cross-sectional  diameter  using  30 
three-dimensional  clay  models  resembling 
the  shapes  of  aggregates  in  our  photographs. 
These  models  were  4-5  cm  long,  radially 
symmetrical,  and  molded  to  approximate 
the  shape  of  the  aggregates  as  seen  from  the 
side  view  in  the  photographs  of  the  natural 
aggregates.  The  projected  area  of  each  mod¬ 
el  was  determined  from  its  maximal  cross- 
sectional  shadow  perpendicular  to  the  polar 
axis.  That  area  was  then  compared  with  the 
calculated  area  of  a  circle  having  a  diameter 
equal  to  the  maximum  cross-sectional  di¬ 
ameter  of  the  model.  The  95%  C.L.  of  the 
mean  on  our  projected  area  measurements 
was  ±  1 9%. 

Volume:  Volume  was  calculated  from  the 
photograph  of  each  aggregate,  assuming  that 
aggregates  were  symmetrical  about  their  po¬ 
lar  axis  and  that  their  volumes  approxi¬ 
mated  ellipsoids  or  spheres.  The  few  irreg¬ 


ularly  shaped  aggregates,  such  as  that  shown 
in  Fig.  le,  were  divided  into  subunits  re¬ 
sembling  spheres,  cylinders,  or  ellipsoids, 
and  total  volume  was  estimated  by  adding 
the  volumes  of  these  subunits.  We  also  es¬ 
timated  the  error  of  calculating  volume  from 
two-dimensional  photographs  using  the  clay 
models.  We  calculated  the  volume  of  these 
clay  models  from  their  two-dimensional 
shadows  projected  in  the  plane  parallel  to 
their  polar  axes  and  compared  these  cal¬ 
culated  volumes  with  the  actual  displace¬ 
ment  volumes  of  the  models  in  water.  These 
trials  yielded  a  95%  C.L.  on  volumes  of 
±30%. 

Porosity;  Porosity  is  the  fraction  of  an 
aggregate  not  occupied  by  solid  matter  and 
was  needed  to  determine  aggregate  density. 
We  calculated  porosity,  P,  directly  from  the 
measurements  of  volume  and  dry  weight 
with  the  equation: 


P  = 


fF/p, 

V 


(5) 


where  W  is  dry  weight  in  g  and  p,  the  density 
of  the  solid  hydrated  matter  within  the  ag¬ 
gregate.  We  assumed  p,  =  1 .23  g  cm“  which 
is  the  wet  density  of  euphausiid  fecal  pellets 
(Komar  et  al.  1981).  Since  zooplankton  fe¬ 
cal  pellets  contain  a  representation  of  many 
of  the  types  of  particles  found  in  marine 
snow  including  diatom  frustules,  clay-min¬ 
eral  particles,  and  intact  cells,  their  wet  den¬ 
sity  seemed  a  reasonable  choice  for  the  den¬ 
sity  of  the  solid  matter  within  aggregates. 
By  comparison,  the  density  of  dry  cellulose 
is  about  1 .5  g  cm~^  and  that  of  living  phy¬ 
toplankton  may  be  equal  to  seawater,  about 
1.025  g  cm"^  (Smayda  1970).  Since  aggre¬ 
gates  were  generally  >  99%  porous,  our  cal¬ 
culated  porosities  were  not  particularly  sen¬ 
sitive  to  the  value  of  p^.  Values  of  p,  0.2  g 
cm“^  higher  or  lower  than  the  value  used 
altered  porosity  of  our  aggregates  only  neg¬ 
ligibly. 


Fig.  1 .  Measurement  of  sinking  rates  of  marine  snow  in  situ.  a.  Dye  ejector  used  by  SCUBA  diver  to  deliver 
neutrally  buoyant  dye  spot  below  aggregate,  b.  Aggregate  formed  from  senescent  diatoms  and  diatom  frustules 
(scale  =  1  cm),  c.  Spherical  mucus  aggregate  formed  from  the  decomposing  house  of  an  appendicularian.  d. 
Comet-shaped  aggregate  of  unknown  origin,  e.  Irregularly  shaped  aggregate  of  unknown  origin  containing  nu¬ 
merous  macrocrustacean  fecal  pellets  (scale  of  c-e  =  I  mm). 
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Excess  density  (Ap);  The  difference  be¬ 
tween  the  density  of  the  aggregate  and  that 
of  the  surrounding  seawater  was  calculated 
from 


Seawater  density  was  determined  from  am¬ 
bient  temperature  and  salinity  according  to 
Mamayev  (1975).  Salinity  was  measured 
with  a  Plessey  Environmental  Systems  lab¬ 
oratory  salinometer  (model  6230N). 

Significant  functional  relationships  be¬ 
tween  the  above  variables  were  determined 
with  methods  of  standard  linear  regression 
(Sokal  and  Rohlf  1969). 

Results 

Characteristics  of  marine  snow— Of  the 
80  aggregates  of  marine  snow  studied  in  the 
field,  68  were  collected  in  the  San  Pedro 
Basin  at  a  salinity  of  33.570%o,  a  tempera¬ 
ture  of  15.0'’C,  and  a  seawater  density  of 
1.02488  g  cm"\  The  remaining  12  were 
from  the  Santa  Barbara  Basin  at  a  salinity 
of  33.642%o,  a  temperature  of  16.0®C,  and 
a  seawater  density  of  1 .02466  g  cm“  Data 
from  the  two  study  sites  were  pooled. 

The  80  aggregates  of  marine  snow  were 
highly  variable  in  size,  shape,  and  general 
appearance.  The  12  aggregates  from  the 
Santa  Barbara  Basin  were  all  flocculent  con¬ 
glomerates  of  living  senescent  diatoms,  par¬ 
ticularly  chain-forming  species,  and  frus- 
tules  (Fig.  lb)  which  formed  following  a 
diatom  bloom  (sf-f- Smetacek  1985).  Aggre¬ 
gates  from  the  San  Pedro  Basin  were  of  di¬ 
verse  origins  and  appearance.  Some  were 
the  remains  of  appendicularian  houses  con¬ 
taining  considerable  gelatinous  mucus  (Fig. 
Ic).  Most  were  aggregates  of  smaller  parti¬ 
cles,  detritus,  and  fecal  pellets.  Shapes  ranged 
from  comets  (Fig.  Id)  to  spheroids  (Fig.  Ic) 
to  oblate  spheres  (Fig.  lb).  Some  of  the  ag¬ 
gregates  had  irregular  shapes  or  contained 
abundant  macrocrustacean  (primarily  eu- 
phausiid)  fecal  pellets  (Fig.  le).  The  shapes 
of  some  aggregates  appeared  to  have  been 
deformed  by  the  flow  of  fluid  around  them 
as  they  sank  (Fig.  lb,  d).  Sinking  aggregates 
maintained  a  stable  orientation  in  the  water. 


and  we  observed  no  twisting  or  rotating  as 
they  sank. 

Diameters  of  the  aggregates  ranged  from 
0.5  to  25.5  mm  and  maximum  lengths  from 
2.4  to  75  mm.  Aggregate  volumes  ranged 
from  3  to  6,000  mm^  and  were  related  to 
approximately  the  cube  of  aggregate  diam¬ 
eter  (Fig.  2A).  Dry  weights  varied  from  2  to 
1 , 1 1 0  jug  aggregate^ '  and  also  increased  sig¬ 
nificantly  with  aggregate  diameter  (Fig.  2B). 
Projected  areas  ranged  from  0. 1  to  5 1 1  mm-. 
Reynolds  numbers  ranged  from  0.4  to  32. 

The  aggregates  studied  had  very  high  po¬ 
rosities  ranging  from  97  to  99.9%  with  a 
mean  of  99.5%.  Porosity  increased  signifi¬ 
cantly  with  increasing  particle  diameter  (Fig. 
2C).  We  plotted  the  log  of  ( 1  —  porosity), 
as  did  Kajihara  (1971),  in  order  to  reveal 
the  large  range  of  porosities  encountered 
(Fig.  2C). 

The  excess  density  (Ap)  of  the  aggregates 
ranged  over  four  orders  of  magnitude  from 
2.2  X  10-Mo  1.3  X  10-5 gem -'(Fig.  2D). 
The  median  was  1.4  x  lO'"*  g  cm-'.  Only 
20%  of  the  aggregates  had  a  Ap  >  1  x  10' 
g  cm”'  (Fig.  2D).  If  we  assume  that  the  in¬ 
terstitial  water  within  the  aggregates  had  the 
same  density  as  the  surrounding  seawater, 
or  1.02488  g  cm"',  the  median  absolute 
density  of  the  aggregates  of  marine  snow 
collected  at  our  study  sites  was  1.02502  g 
cm"'.  Excess  density  was  highly  correlated 
with  aggregate  size  (Fig.  2D)  despite  the 
variable  nature  of  the  primary  particles 
comprising  the  aggregates— including  vary¬ 
ing  proportions  of  fecal  pellets,  mucus, 
molts,  frustules,  and  other  components.  Al¬ 
though  it  appeared  that  diatom  floes  might 
have  a  different  excess  density  to  size  rela¬ 
tionship  than  other  types  of  marine  snow, 
the  slope  of  the  regression  line  describing 
that  relationship  was  not  significantly  dif¬ 
ferent  (ANCOVA  test  for  equality  of  slope: 
Sokal  and  Rohlf  1 969)  from  the  slope  of  the 
regression  line  describing  the  excess  density 
to  size  relationship  for  all  of  the  other  ag¬ 
gregates  (Fig.  2D).  The  data  set  for  diatom 
snow  was  small  with  high  variability,  how¬ 
ever,  so  additional  data  may  reveal  a  sig¬ 
nificant  difference. 

Settling  velocities— T'xguTt  3  displays  the 
settling  velocities  of  80  aggregates  of  marine 
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DIAMETER  (mm) 


Fig.  2.  Size  characteristics  of  marine  snow.  A.  Aggregate  volume  as  a  function  of  diameter  (regression 
coefficient,  rc,  =  0.97;  P  <  0.001).  B.  Aggregate  dry  weight  as  a  function  of  diameter  (rc  =  0.71;  Z’  <  0.001).  C. 
Aggregate  porosity  as  a  function  of  diameter;  O— diatom  floes;  •—all  other  marine  snow  regardless  of  origin 
(rc  =  0.79;  P  <  0.001).  D.  Aggregate  excess  density  as  a  function  of  diameter  (rc  =  0.79;  P  <  0.001),  symbols 
as  in  panel  C. 


snow  in  situ  as  a  function  of  various  particle 
characteristics.  Mean  settling  rate  was 
74±39  m  d“'.  Sinking  speed  increased  ex¬ 
ponentially  with  particle  diameter  (Fig.  3A). 
Settling  velocities  increased  with  the  in¬ 
creasing  ratio  of  volume  to  projected  area, 
as  predicted  by  settling  theory  (Fig.  3B). 

Settling  rate  in  situ  increased  exponen¬ 
tially  with  aggregate  dry  weight  (Fig.  3C). 
We  compared  the  size-specific  settling  rates 
of  marine  snow  determined  in  the  labora¬ 
tory  with  rates  for  similarly  sized  and  shaped 
aggregates  determined  in  situ  with  dry  weight 
as  an  accurate  measure  of  aggregate  size. 
Despite  minimal  handling  of  particles,  set¬ 
tling  velocities  measured  in  the  laboratory 
were  consistently  higher,  by  up  to  four  times, 
than  those  of  similarly  sized  aggregates 
measured  in  situ  (Fig.  3C).  No  significant 


statistical  correlation  could  be  found  be¬ 
tween  sinking  rate  and  dry  weight  of  aggre¬ 
gates  studied  in  the  laboratory. 

Although  settling  theory  predicts  that  the 
sinking  rate  of  an  object  settling  in  a  fluid 
is  a  function  of  the  excess  density  of  the 
object,  our  data  did  rot  yield  a  significant 
relationship  between  excess  aggregate  den¬ 
sity,  as  calculated  by  our  methods,  and  sink¬ 
ing  rate  (Fig.  3D). 

We  measured  or  derived  all  of  the  di¬ 
mensional  terms  in  the  force  balance  Eq.  1 
for  a  settling  object.  Thus,  we  can  directly 
calculate  the  drag  coefficient,  C„,  for  each 
sinking  aggregate.  Calculations  of  C/,  for  our 
nonspherical  aggregates  can  provide  insight 
into  the  effects  of  shape  and  other  variables 
that  potentially  alter  the  settling  behavior 
of  marine  snow  relative  to  sinking  spheres. 


346 


Alldredge  and  Gotschalk 


DIAMETER  (mm) 


VOLUME/CROSS-SECTIONAL  AREA  (cm) 


10°  1o’  10^  10° 


DRY  WEIGHT  (pg) 


EXCESS  DENSITY  (g  cm"'’) 


Fig.  3.  Sinking  rates  of  marine  snow  in  situ  (•)  as  a  function  of:  A— aggregate  diameter  (regression  coefficient, 
rc.  =  0.61:  P  <  0.001);  B— aggregate  volume  divided  by  projected  area  (rc  =  0.42;  P  <  0.001);  C— aggregate 
dry  weight  (rc  =  0.60;  P  <  0.001)  (O— sinking  velocities  determined  in  the  laboratory);  D— aggregate  excess 
density. 


Figure  4A  presents  as  a  function  of 
Reynolds  number.  Also  plotted  on  Fig.  4A 
are  the  predicted  curves  for  spheres  assum¬ 
ing  Stokes  settling  (C„  =  24/Re)  and  an  em¬ 
pirically  determined  relationship. 


24 

-  + 

Re 


6 

1  -F  (Re)"’ 


+  0.4, 


(7) 


for  solid  spheres  at  higher  Reynolds  num¬ 
bers  (White  1974). 

The  drag  coefficients  of  marine  snow  were 
consistently  higher  than  those  of  spheres  of 
equivalent  Reynolds  number  (Fig.  4A),  in¬ 
dicating  that  marine  snow  sank  more  slowly 
than  equivalent  spheres.  Spheres  sink  more 
rapidly  than  most  other  shapes  (Simpson 
1982).  Moreover,  for  nonspherical  parti¬ 
cles,  C„  is  a  function  of  both  Re  and  shape. 
Therefore,  we  investigated  shape  as  a  factor 
affecting  settling  velocity. 


We  hypothesized  that  the  coefficient  of 
form  resistance,  k,  would  also  increase  with 
increasing  nonsphericity  for  marine  snow 
particles.  We  determined  the  Corey  shape 
factor,  CSF,  assuming  that  the  particles  were 
symmetrical  about  their  polar  axis;  thus  the 
two  smallest  axial  diameters,  D,  and  £>„  were 
equal  for  prolate  aggregates  while  D,  and  D/ 
were  equal  for  oblate  aggregates.  We  cal¬ 
culated  t/„  the  settling  velocity  of  a  sphere 
of  equivalent  volume  and  density  to  each 
aggregate,  with  the  force  balance  Eq.  2.  In 
solving  this  equation  for  equivalent  spheres, 
we  used  the  equation  for  C^,  empirically  de¬ 
rived  by  White  (Eq.  7)  for  spheres  outside 
the  Stokes  range.  Since  excess  density  was 
size-specific,  we  used  the  empirically  de¬ 
rived  function  of  Fig.  2D  to  calculate  the 
excess  density  of  sinking  spheres  as  a  func¬ 
tion  of  sphere  diameter.  Substitution  oi  these 
parameters  into  the  force  balance  Eq.  2  en- 
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abled  us  to  obtain  a  solution  for  the  settling 
velocity  of  a  sphere  with  a  volume  and  den¬ 
sity  equivalent  to  each  aggregate  of  marine 
snow. 

We  found  no  significant  relationship  be¬ 
tween  k  and  the  shapes  of  the  aggregates 
(Fig.  4B).  We  also  found  no  relationship 
when  the  analysis  was  restricted  to  oblate 
aggregates  or  prolate  aggregates  only.  Ag¬ 
gregates  whose  shapes  varied  greatly  from 
spherical,  including  long  comets,  had  values 
of  k  very  similar  to  those  of  nearly  spherical 
particles.  Some  near-spherical  aggregates 
sank  considerably  slower  than  equivalent 
spheres,  with  values  of  /c  S  3,  rather  than 
the  predicted  value  of  1.  Previous  studies 
on  the  effects  of  shape  on  predict  that 
particles  with  highly  nonsperhical  shapes 
would  have  higher  drag  coefficients  (Komar 
and  Reimers  1978),  but  we  found  no  sig¬ 
nificant  relationship  between  shape  and  Q, 
(Fig.  4C). 

The  C[)  of  marine  snow  was  higher  than 
predicted  for  sinking  spheres  of  equivalent 
Reynolds  number,  since  the  settling  velocity 
of  marine  snow  was  slower  than  that  of 
equivalent  spherical  particles.  Although  we 
found  no  relationship  between  either  drag 
coefficient  or  the  coefficient  of  form  resis¬ 
tance  and  the  Corey  shape  factor,  shape 
might  still  partially  explain  the  high  Co  ob¬ 
served.  Other  factors  including  permeabil¬ 
ity  and  surface  roughness  might  also  in¬ 
crease  drag. 

In  order  to  elucidate  these  issues,  we  plot¬ 
ted  the  sinking  velocity  of  particles  of  sev¬ 
eral  different  shapes  and  porosities  having 
volumes  and  excess  densities  equivalent  to 
those  of  our  aggregates  (Fig.  5).  The  equa¬ 
tions  for  the  settling  velocities  of  a  sphere, 
and  a  prolate  ellipsoid  with  an  aspect  ratio 
of  4 :  1  [prolate  ellipsoid  (a)  based  on  Stokes’ 
law,  Fig.  5]  were  obtained  from  Lerman  et 
al.  (1974).  The  settling  velocities  of  equiv¬ 
alent  spheres  at  Reynolds  numbers  outside 
the  Stokes  range  (Re  >  0.5)  were  obtained 
from  White  (1974)  as  described  previously. 
The  empirically  derived  curve  for  the  set¬ 
tling  velocity  of  equivalent  prolate  ellip¬ 
soids  with  aspect  ratio  of  4 : 1  at  Re  >  0.5 
[prolate  eUipsoid  (b).  Fig.  5]  was  obtained 
by  solving  Eq.  4  for  the  settling  velocity,  V, 
of  a  prolate  ellipsoid  by  using  the  settling 
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Fig.  4.  Drag  coefficient  (C„)  and  effects  of  shape 
(O— diatom  snow;  •—marine  snow  of  all  other  origins). 
A.  Co  vs.  Reynolds  number  (Re);  White— Co  of  sinking 
spheres  determined  empirically  from  White  (1974); 
Stokes— Co,  as  predicted  by  Stokes’  equation.  B.  The 
dynamic  shape  factor,  k,  as  a  function  of  the  Corey 
shape  factor  (CSF).  CSF  should  equal  1  for  perfect 
spheres  (see  text  for  definitions).  C.  Drag  coefficient  as 
a  function  of  the  Corey  shape  factor. 
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Fig.  5.  Settling  velocities  of  marine  snow  (•)  and 
of  various  spheres  and  prolate  ellipsoids  with  volumes 
and  excess  densities  equal  to  those  of  the  aggregates 
studied.  Dashed  lines  are  theoretical  settling  rates  pre¬ 
dicted  by  Stokes’  (“Stokes  sphere”)  or  variations  of 
Stokes’  law  for  prolate  ellipsoids,  prolate  4: 1  (a).  Solid 
lines  are  settling  velocities  derived  from  empirical 
.studies.  "White  sphere,”  solid  spheres  according  to 
White  (1974).  Porous  spheres  predicted  from  empirical 
data  of  Matsumoto  and  Suganuma  (1977).  Prolate  4;  1 
(b)  predicted  from  empirical  data  of  White  ( 1 974)  with 
k  values  from  McNewn  and  Malaika  (1950)  to  correct 
for  nonspherical  shape. 


velocities  of  sinking  spheres  (U,)  from  White 
(1974)  and  k  values  for  a  prolate  ellipsoid 
of  4 :  1  aspect  latio  from  McNown  and  Ma¬ 
laika  (1950)  who  empirically  determined  k 
up  to  Re  of  ~10. 

In  order  to  estimate  the  potential  effects 
of  permeability  on  settling  velocities  of  mu- 
rine  snow,  we  used  data  on  permeable  floes 
of  steel  wool  ranging  in  porosity  from  90  to 
99.7%  and  in  size  from  3  to  100  mm  studied 
by  Matsumoto  and  Suganuma  (1977).  The 
settling  velocities  of  permeable  steel  wool 
floes  having  volumes  and  densities  equiv¬ 
alent  to  those  of  the  marine  snow  were  cal¬ 
culated  by  correcting  sinking  rates  of  solid 
equivalent  spheres  as  determined  with 
White’s  equation,  using  the  size-specific 
correction  factor,  F,  developed  empirically 
by  Matsumoto  and  Suganuma  (1977). 

Figure  5  demonstrates  that  aggregates  sink 
considerably  more  slowly  than  do  spheres 
and  4 ;  1  prolate  spheres  of  equivalent  vol¬ 
ume  and  excess  density  as  predicted  with 
Stokes’  law.  Aggregates  sink,  on  average, 
only  about  half  as  fast,  however,  as  would 
be  expected  for  equivalent  spheres  at  Reyn¬ 
olds  numbers  >  1 .  Porosity  appears  to  in¬ 


crease,  rather  than  decrease,  sinking  veloc¬ 
ities  within  the  range  studied.  The  porous 
spheres  of  Matsumoto  and  Suganuma  (1977) 
had  rough,  convoluted  surfaces  analogous 
to  the  surfaces  of  marine  snow  aggregates. 
Such  convoluted  surfaces  might  increase 
skin  friction  and  thus  total  drag  on  the 
particle.  Neither  surface  roughness  nor  per¬ 
meability,  however,  appears  to  have  in¬ 
creased  the  drag  of  Matsumoto  and  Suga- 
numa’s  particles.  Masliyah  and  Polikar 
(1980)  also  investigated  the  settling  veloc¬ 
ities  of  permeable  spheres  of  milled  plastic 
up  to  2.5  cm  in  diameter  and  at  Re  up  to 
120.  These  spheres  were  97%  porous  with 
highly  convoluted  and  prickly  surfaces. 
Their  settling  velocities  also  deviated  only 
slightly  from  the  settling  velocities  of  equiv¬ 
alent  impermeable  spheres.  Thus  perme¬ 
ability  and  surface  roughness  probably  have 
relatively  little  effect  on  the  settling  veloc¬ 
ities  of  marine  snow.  Logan  and  Hunt  (1987) 
concluded  that  surface  roughness  and  flow 
through  microbial  aggregates  probably  did 
not  contribute  to  observed  deviations  in  set¬ 
tling  velocity  compared  to  predictions  for 
impermeable  spheres. 

Discussion 

Many  different  methods  have  been  used 
to  estimate  the  sinking  rates  of  marine  snow. 
Our  sinking  rates  tended  to  be  lower  than 
those  estimated  from  time-lapse  photogra¬ 
phy  of  the  seafloor  (100-150  m  d  Billett 
et  al.  1983;  Lampitt  1985)  and  lower  than 
those  measured  by  ourselves  and  others  in 
settling  chambers  in  the  laboratory  (Silver 
and  Alldredge  1981;  Taguchi  1982;  Gorsky 
et  al.  1983).  Collection  and  gentle  transfer 
of  fragile  aggregates  in  the  laboratory  ap¬ 
parently  result  in  slight  collapse  of  the  ag¬ 
gregate,  with  concurrent  decreases  in  po¬ 
rosity,  increases  in  density,  and  increases  in 
settling  velocity.  Our  sinking  rates  were 
similar  to  the  range  of 43-95  m  d  '  reported 
by  Shanks  and  Trent  (1980),  who  measured 
sinking  rates  of  aggregates  in  large  cylinders 
in  situ.  In  general,  we  observed  similar  sink¬ 
ing  rates  to  those  calculated  by  Asper  ( 1 986) 
from  sediment  trap  data  and  particle  abun¬ 
dances  in  the  water  column,  although  we 
never  observed  rates  as  low  as  the  1  m  d  ' 
he  reported  from  deep  water  of  the  Panama 
Basin.  The  highly  porous,  flocculent  nature 
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of  the  material  he  observed  may  explain  its 
slow  sinking  rate. 

Previous  investigations  of  the  settling  ve¬ 
locities  of  natural  aggregated  matter  as  a 
function  of  aggregate  diameter  vary  some¬ 
what  from  our  results.  Kajihara  (1971)  found 
U=  160rf"^^and Gibbs(1985) reported  V  = 
respectively  for  reaggregated  and 
estuarine  floes  of  <  1  mm;  we  report  U  = 
for  marine  snow  aggregates  of  larger 

size. 

Previous  investigations  of  the  porosity  and 
excess  density  of  natural  and  manmade  floes 
have  yielded  results  similar  to  those  re¬ 
ported  here.  Assuming  Stokes’  law,  Kaji¬ 
hara  (1971)  claculated  the  porosity  of  nat¬ 
ural  marine  floes  <  1  mm  from  settling 
velocities  measured  in  the  laboratory.  He 
found  that  porosity  increased  with  increas¬ 
ing  particle  size,  although  he  reported  po¬ 
rosities  lower  than  ours  for  similarly  sized 
particles.  His  particles  were  formed  by  reag¬ 
gregation  of  smaller  component  particles  and 
were  probably  more  dense  than  natural  floes. 
Gibbs  (1985)  found  that  excess  density  of 
estuarine  floes  <  1  mm  decreased  with  in¬ 
creasing  floe  size.  Tambo  and  Watanabe 
(1979)  also  found  decreasing  excess  density 
with  increasing  size  for  clay-aluminum  floes 
in  the  laboratory,  although  the  excess  den¬ 
sity  of  these  inorganic  floes  was  about  three 
times  that  of  marine  snow  aggregates  of 
equivalent  size.  Our  empirical  result  that 
porosity  of  marine  snow  increases  in  pro¬ 
portion  to  d'*  corroborates  the  indepen¬ 
dently  derived  assumption  of  Logan  and 
Hunt  (1987)  that  the  porosity  of  microbial 
floes  increases  in  proportion  to  d  ' 

We  observed  a  large  range  and  high  scat¬ 
ter  of  settling  velocities  when  plotted  against 
diameter  and  dry  weight  — parameters  that 
could  be  measured  with  minimal  error.  This 
scatter  most  likely  arose  from  error  in  mea¬ 
suring  sinking  rates  in  situ  (variations  in 
exact  sinking  distance,  etc.)  and  from  the 
highly  heterogeneous  nature  of  the  aggre¬ 
gates  themselves.  Aggregate  size,  shape, 
composition  (i.e.  presence  of  fecal  pellets, 
mucus,  diatoms,  etc.),  porosity,  and  rigidity 
were  all  highly  variable  in  the  population 
of  natural  and  randomly  chosen  particles 
studied  here,  resulting  in  a  large  range  of 
sinking  rates.  Plots  of  settling  velocity  vs. 
size  for  natural  microscopic  particles  of  sed¬ 


iment  of  uniform  composition  show  scatter 
similar  to  that  which  we  observed  for  mac¬ 
roscopic  aggregates  (Hawley  1982). 

Nonspherical  shapes  result  in  decreased 
settling  velocities  of  aggregates  in  all  the  cases 
examined  in  Fig.  5.  In  fact,  when  nonspher¬ 
ical  shape  is  taken  into  account,  macro¬ 
scopic  aggregates  appear  to  sink  within  the 
range  predicted  for  similar  nonspherical 
particles  sinking  outside  the  Stokes  range 
[Fig.  5,  curve  4: 1  prolate  (b)].  Thus,  devia¬ 
tions  in  the  settling  rate  of  marine  snow 
from  that  of  spheres  at  Reynolds  numbers 
>  I  may  result  partly  from  deviations  from 
nonspherical  shape.  Why  then  did  we  not 
obtain  significant  relationships  between  set¬ 
tling  velocity  and  either  shape  or  excess  den¬ 
sity? 

The  coefficient  of  form  resistence  {k)  and 
settling  velocity  both  depend  on  excess  den¬ 
sity.  Excess  density  is  directly  dependent  on 
the  primary  particle  density  (see  Eq.  6).  We 
assumed  that  this  primary  particle  density 
was  a  constant  for  all  aggregates  and  equal 
to  the  wet  density  of  euphausiid  fecal  pel- 
lets—the  only  density  value  we  could  find 
in  the  literature  for  a  natural  marine  particle 
with  components  similar  to  those  of  marine 
snow.  This  assumption  probably  resulted  in 
sufficient  error  to  obscure  the  expected  re¬ 
lationship  between  sinking  and  shape  and 
excess  density  of  marine  snow. 

Assumption  of  a  constant  primary  par¬ 
ticle  density  for  all  marine  snow  aggregates 
is  problematic  for  two  major  reasons.  First, 
the  primary  panicles  composing  marine 
snow  aggregates  are  not  homogeneous.  Dia¬ 
tom  frustules,  clay  mineral  particles,  fecal 
pellets,  living  phytoplankton,  protozoa, 
crustacean  carapaces,  and  the  various  un¬ 
identifiable  debris  making  up  marine  snow 
each  Imvc  their  own  densities.  Each  aggre¬ 
gate  is  composed  of  varying  percentages  of 
these  and  other  components  yielding  its  own, 
unique  total  excess  density.  To  accurately 
determine  the  excess  density  of  each  aggre¬ 
gate  we  would  have  to  know  not  only  the 
exact  composition  and  size  of  each  com¬ 
ponent  particle  comprising  it,  but  also  the 
densities  of  those  components.  Reliable  data 
of  this  sort  are  not  currently  available  for 
any  of  the  complex  types  of  natural  particles 
composing  marine  snow. 

Second,  the  problem  of  determining  the 
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density  of  the  component  particles  in  ma¬ 
rine  snow  is  further  confounded  by  the  siz¬ 
able  percentage  of  living  organisms  com¬ 
prising  some  aggregates  (Silver  et  al.  1978). 
Many  of  these  organisms  can  regulate  their 
individual  buoyancies.  Diatom  floes  rep¬ 
resent  an  extreme  type  of  marine  snow  where 
>90%  of  the  particle  may  be  living  cells 
(Alldredge  unpubl.).  The  density  of  diatoms 
varies  with  their  physiological  state  and  even 
with  the  time  of  day  (Eppley  et  al.  1967; 
Smayda  1970).  When  nutrient  replete,  dia¬ 
tom  cells  may  be  neutrally  buoyant.  When 
nutrient  stressed  their  excess  density  ranges 
up  to  0.08  g  cm”^  (Eppley  et  al.  1967).  It  is 
probable  that  we  overestimated  the  excess 
density  of  the  diatom  floes  because  we  as¬ 
sumed  they  had  a  constant  primary  particle 
density  similar  to  other  aggregates  that  con¬ 
tain  a  much  lower  proportion  of  living  cells. 
This  assumption  may  have  produced  the 
apparent  deviations  in  porosity,  excess  den¬ 
sity,  and  drag  coefficient  observed  for  this 
type  relpfive  to  other  types  of  marine  snow. 
.Aggregates  are  not  composed  solely  of  in¬ 
animate  particles,  and  thus,  we  would  ex¬ 
pect  their  overall  density  to  depend,  in  part, 
on  the  density  of  their  living  components. 
This  density  may  vary  on  a  scale  as  short 
as  hours,  depending  upon  the  nature  of  the 
community  on  the  snow  and  its  physiolog¬ 
ical  state.  (Tlearly  additional  research  on  the 
density  of  natural  marine  particles  is  needed 
to  refine  the  relationships  between  sinking 
and  aggregate  characteristics  reported  here. 

Relationships  between  sinking  and  both 
excess  density  and  shape  may  also  have  been 
obscured  by  high  scatter  in  the  data  and  by 
the  compounding  of  error  in  the  measure¬ 
ments.  Excess  density  is  dependent  on  vol¬ 
ume  and  dry  weight— each  with  its  own 
corresponding  errors.  The  odd  finding  that 
some  of  the  values  for  the  drag  coefficients 
of  marine  snow  are  lower  than  those  pre¬ 
dicted  by  Stokes’  law  may  also  arise  from 
compounding  of  random  error  since  the  cal¬ 
culation  of  drag  coefficient  depends  on  sev¬ 
eral  empirically  determined  variables  in¬ 
cluding  settling  velocity,  projected  area,  and 
excess  density. 

Data  provided  herein  may  allow  esti¬ 
mation  of  the  characteristics  and  flux  of  nat¬ 
ural  particles  based  on  their  size  distribu¬ 
tions  in  situ  generated  from  photographs 


and  survey  cameras  (see  Asper  1987;  John¬ 
son  and  Wangersky  1 985).  Estimates  of  par¬ 
ticulate  flux  based  on  collection  of  particles 
in  the  water  column  with  in  situ  pumps  have 
already  provided  considerable  information 
about  flux  (Bishop  et  al.  1977,  1980).  Care 
must  be  taken  in  calculating  flux  from  pre¬ 
dicted  settling  rates  and  measured  particle 
abundances  in  the  water  column,  however. 
Just  because  large  particles  sink  does  not 
necessarily  mean  that  they  sink  out  of  the 
water  column  at  the  predicted  rate.  Accu¬ 
mulation  of  marine  snow  has  been  docu¬ 
mented  at  a  well-developed  pycnocline  in 
the  subtropical  Atlantic  (Alldredge  and 
Youngbluth  1985).  Marine  snow  accumu¬ 
lating  at  the  top  of  a  density  discontinuity 
could  be  reinjected  back  into  the  mixed  lay¬ 
er  by  wind-mixing  events.  Reinjection  of 
some  proportion  of  the  marine  snow  pop¬ 
ulation  back  into  the  mixed  layer  would  re¬ 
sult  in  accumulation  of  marine  snow  in  sur¬ 
face  waters  and  potential  residence  times  for 
these  particles  longer  than  would  be  pre¬ 
dicted  based  on  their  sinking  rates  alone. 
Macrocrustacean  fecal  pellets,  particles  with 
settling  velocities  and  excess  densities  con¬ 
siderably  greater  than  those  of  marine  snow, 
are  known  to  have  residence  times  in  sur¬ 
face  waters  up  to  10  d  or  more  (Alldredge 
etal.  1987).  Moreover,  lateral  advection  ap¬ 
pears  responsible  for  high  snow  abundances 
at  certain  depths  in  the  Panama  Basin  (As¬ 
per  1986);  aggregates  are  carried  horizon¬ 
tally  as  well  as  sinking  vertically.  Since  flux 
calculations  often  assume  that  the  particles 
immediately  sink  out  of  the  system,  calcu¬ 
lations  of  particulate  flux  based  on  particle 
abundances  and  size  distributions  in  situ 
and  on  predicted  settling  velocities  could 
overestimate  flux  to  the  ocean  interior  and 
to  the  seafloor. 
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